Is the island universe model consistent with observations? 
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We study the island universe model, in which initially the universe is in a cosmological constant 
sea, then the local quantum fluctuations violating the null energy condition create the islands of 
matter, some of which might corresponds to our observable universe. We examine the possibility 
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that the island universe model is regarded as an alternative scenario of the origin of observable 
O ' universe. 
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Recent observations of the cosmic microwave background (CMB) have implied that the 
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inflation [jj, |2J is a very consistent cosmological scenario. The inflation stage took place at 
the earlier moments of the universe, which superluminally stretched a tiny patch to become 
our observable universe today, and ended with a period of reheating that bring the universe 
back to the usual FRW evolution. The quantum fluctuations in the inflation period are drew 
out the horizon and, when reentering into the horizon during radiation/matter domination 
after the end of inflation, are responsible for the formation of cosmological structure. 

However, the inflation has still its own problem of fine tuning, in addition it is not also 
able to address the initial singularity problem, the cosmological constant problem, super 
Planck fluctuations problem, as was described in Ref. (^J. Thus it seems desirable to 
search for an alternative to the inflation model. Recently, Dutta and Vachaspatian have 
proposed a model ^], in which initially the universe is in a cosmological constant ( A) sea, 
then the local quantum fluctuations violating the null energy condition (NEC) 0, Q create 
the islands of matter, some of which might corresponds to our observable universe. The 
headstream of this model in some sense may be backward to the eternal inflation [3, 0, 
and especially the recycling universe proposed by Garriga and Vilenkin j|, in which the 
thermalised regions asymptotically approaching dS spacetime may be fluctuated and recycled 
back to the false vacuum and then the nucleated false vacuum region will serve as a seed 
for a new (eternally) inflating domain. The difference of island universe is that after the 
NEC violating fluctuation is over, instead of the inflation followed, the thermalisation occurs 
instantly and the radiation fills the volume rapidly. Though the island universe model brings 
us an interesting environment, people are still able to doubt whether the island universe 
model is consistent with our observations. 

Whether we can live in "Islands"? We firstly attempt to phenomenally (semiclassically) 
understand it, following closely Ref. Q|. We set Sir/m 2 = 1 and work with the parameter 
e = —h/h 2 , where h = a/a is the local Hubble constant. The "local" means here that the 
quantities, such as a and h, only character the value of the NEC violating region. In the 
island universe model, the scale of the NEC violating region is required to be larger than 
that of present horizon, which has been clearly stated by Dutta and Vachaspatian |4j, and 
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for relevant discussions. This 



result provides the initial value of local evolution of a and h. The e can be rewritten as 
e ~ tkil! 1 ' ^hus in some sense e actually describes the change of h in unit of Hubble time. 
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During the NEC violating fluctuation, h > 0, thus e < can be deduced. The time of this 
fluctuation can be given by 

f , r dh 1 

where the subscript i and e denote the initial and end value of the NEC violating fluctuation 
respectively. In this model h\ = h\ ~ A, where the subscript denotes the present value. 
The reasonable and simplest selection for the local evolution of scale factor a(t), due to 
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h > 0, is |10L Il7| a(t) ~ (— t) n , where i is from — oo to 0_, and n is a negative constant. 
We have h = n/t, and thus e = 1/n. To make T — > in which the NEC violating fluctuate 
can be so strong as to be able to create the islands of our observable universe from 
|e| — > oo is required, which results in n — > 0_. Thus though the change of h is large, the 
expanding proportion of the scale factor is very small. Further in some sense before and after 
the fluctuation a is nearly unchanged, which can be also seen from following discussions. 
In the conformal time dr\ = dt/a, i.e. —7] ~ (— t)~ n+1 , we obtain 



a(ij) ~ (—77) 1-™ = (-77) -1 (2) 

h = ^= (e-l)a V ' (3) 
where the prime denotes the derivative with respect to 77. The perturbations leaving the hori- 
zon during the NEC violating fluctuation might reenter the horizon during radiation/matter 
domination after the end of fluctuation, which are responsible for the structure formation of 
our observable universe. The efolding number which measures the quantity of perturbation 
leaving the horizon during the NEC violating fluctuation can be defined as 

= (4) 
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From (J2J and Q, we obtain 



Thus for the constant e, we have 
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— = (— = e i-«. (6) 

CLi Ct e ll e 

We can see that for the negative enough e, the change Aa/a = (a e — a^jai ~ N / (1 — e) of 
a will be very small. 
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FIG. 1: The sketch of evolution of local In (1/ ah) with respect to the scale factor In a, in which 
the yellow region denotes the NEC violating fluctuation and the red lines denote the evolution of 
primordial perturbation modes. The perturbation modes exit the Hubble horizon during the NEC 
violating fluctuation and then reenter the horizon during the radiation/matter domination at late 
time. 



FIG. 2: The Penrose diagram of creation and evolution of island universe, in which the blue line 
is the event horizon of A sea, the green line and red lines denote the evolution of Hubble scale and 
primordial perturbation modes respectively. The region between black dashed lines denotes the 
NEC violating fluctuation. The dashed line of t = is the thermalisation surface, after which the 
NEC violating fluctuation is over, the local island universe starts the FRW evolution, and with the 
lapsing of time finally will return to the blue A sea again. 

Taking the logarithm in both sides of (f5|). we obtain 



We plot Fig.l to further illustrate the island universe model, in which In (1/ah) is regarded 
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as the function of In a. In the NEC violating region, n — > 0_, thus the slop (1 — n)/n — > — oo, 
which makes the evolution of In (1/ah) during the fluctuation correspond to a nearly vertical 
line in Fig.l. The (JJJ) can be also shown and actually applied for the expansions with 
arbitrary constant n. Thus during the A, radiation and matter domination, we have the 
slope —1, 1 and 1/2 respectively, which have been reflected in the Fig.l. 

The NEC violating fluctuations in the A sea with the observable value of cosmological 
constant create some thermalised matter islands, which subsequently evolute as the usual 
FRW universe. The radiation and matter will be diluted with the expansion of island and 
eventually this part of volume will return to the A sea again. The total evolution may be 
depicted in the Penrose diagram of Fig. 2. We can see that there are always some perturba- 
tions (b) never reentering and remaining outside the Hubble scale after their leaving from 
it during the NEC violating fluctuation, which means that a part of island is permanently 
inaccessible to any given observer in the island. This may be reanalysed as follows. The 
wavelength of perturbations grows with the scale factor ~ a, thus for the perturbations with 
the present horizon scale, we have ao^o = ctch c , where the subscript c denotes the value at 
the time when this perturbation leaves the horizon during the NEC violating fluctuation. 
We define 

as the efolding number when the perturbation with the present horizon scale leaves the 
horizon during the NEC violating fluctuation, and from (j3|, we have 

Af ei -K c = \nA. (9) 

di Hi 

Thus for the island universe model, in which hi = h and a ^> a,i > l//tj, we obtain 
M e i > J\f ec . This implies that as long as the islands created are suitable for our existence, 
the efolding number required to solve the horizon problem of FRW universe may be always 
enough, independent of the energy scale of thermalisation. 

Then we discuss the the primordial perturbations from the island universe model. Dutta 
and Vachaspatian j^J have shown that the spectrum from the quantum fields other than the 
NEC violating field can be scale invariant, but the amplitude of perturbation is too small 
to seed the structure of observable universe. However, they in their calculations regarded 
the NEC violating region as an instant link between the A sea and FRW evolution, which 
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to some extent leads that the effect of NEC violating region on the perturbation spectrums 
is lose. Thus we here will focus on the NEC violating region. 

To simplification, we adopt the work hypothesis in which the NEC violating fluctuation 
behaves like the phantom energy. The simplest implementation of phantom energy is a 
scalar field with the reverse sign in its dynamical term, in which the NEC p + p = —ip 2 < 
is violated. To make p > 0, the potential V(<p) > <p 2 /2 of phantom field is also required. 
Thus in this case the local evolution of h can be written as h 2 ~ — oi 2 /2 + V((p). This in 

n 

some sense is equal to introducing a creation field, see Ref. [18]. Thus we have 

ii P + P 2 nm 

w + 1 = = 7T, . = -e — > — oo. (10) 

p -p 2 /2 + V(p) 3 V ; 

Due to V(<p) > <p 2 /2, this inequality actually requires the denominator ~ 0, which gives 
V((p) ~ 2 /2. Further, combining <p 2 /2 = h = \e\h 2 , we may approximately deduce 

ip 2 /2 ~ V(<p) ~ \e\h 2 . (11) 

This operation dose not mean that the phantom field is actually required in the island 
universe model. The aim that we appeal to the phantom field here is only to phenomenally 
simulate the NEC violating fluctuation, which may be convenient to calculate the primordial 
perturbation spectra from the NEC violating region. But it maybe also a possibility that the 
phantom field only exists instantly for producing the NEC violating fluctuation and after 
the end of fluctuation, the phantom field thermalises into the radiation, see Ref. Q] for a 
different scenario. 

In the momentum space, the equation of motion of gauge invariant variable Uk, which is 
related to the Bardeen potential $ by Uk = a^^/tp', is 1 

< + (k 2 ~ — )u k = 0. (14) 



1 In the momentum space, the equation of motion of gauge invariant variable Vk, which is related to the 
tensor perturbation by Vk = ahk, is 



a 



v k =0. (12) 



We have, from a" /a = (2 — e)/(e— l) 2 ij 2 . In the regime ki] — > oo, the mode Vk ~ ^k) 1 ' 2 can ^ e taken 
as the initial condition. In long wave limit, the expansion of the Bessel functions to the leading term of 
k gives 

k 3 / 2 v k /a ~ fc-'+ 3 / 2 (13) 



where v = y/{2 - e)/(e - l) 2 + 1/4. 
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We have, from (J2J), 

t - (V*)" - e mi 

2 (l/a) (e - 1) V 1 ' 

The general solutions of this equation are the Bessels functions. In the regime krj — > oo, the 
mode Uk are very deep in the horizon of A, see Fig.l and Fig.2. Thus Uk ~ ^k)^/ 2 can be 
taken as the initial condition. In the regime krj — ► 0, the mode are far out the horizon, 
and become unstable and grows. In long wave limit, the expansion of the Bessel functions 
to the leading term of k gives 



k 3/2 u k ~ 93/2 . ( ^ rn — -A-^r +l/2 , (i6) 

2 6 / 2 sm [ttv)i (1 — v) 2 

where v = \/e/(e — l) 2 + 1/4. We can see that if v ~ 1/2, the spectrum of will be nearly 
scale invariant. This requires 

Thus for the NEC violating fluctuation with e — > — oo 2 , we will have the nearly scale 



invariant spectrum 



id- 



32j $ = C + hV/a, where 



In krj — > 0, the Bardeen potential can also be rewritten as 
C and T> are constants dependent on the mode k. In the inflation model in which e ~ 0, 
h is nearly unchanged while a increases exponentially, thus T> is a decay mode, and $ is 
dominated by the constant mode C, and thus is nearly constant in superhorizon scale. This 
is the reason why people can briefly obtain the amplitude of perturbation reentering into 
the horizon by calculating the amplitude of perturbation leaving the horizon during the 
inflation. But for e — > — oo, the case is just inverse, a is nearly unchanged while h increases 
rapidly, thus instead of being regarded as the decay mode, T> is a growing mode, which can 
be seen from (|18|). and compared with the constant mode C, it will dominate the Bardeen 
potential <3>, which means that the amplitude of perturbation in the superhorizon scale will 
be still increasing during the NEC violating fluctuation, while during the radiation/matter 
domination after the thermalisation, T> will not increase any more but decrease, and thus the 



2 We notice from <|1 7f) that for e ~ 0_ in which the scale factor expands exponentially, the spectrum is 
also nearly scale invariant. This is the so called phantom inflation proposed and discussed in detail in 
Ref. 0, see also Refs. 0, li^, I23I Q • This result is actually also a direct reflection of the dualities of 
primordial perturbation spectra, in which the spectrum index of $ is invariant under the change e — ■> 1/e, 
which has been pointed out in Refs. [3)E3| and further studied in Refs. 0,E3jE3j an d see 8 ^ so Rcfs- 
29, 30, 3 1 1 for other discussions on the dualities. 
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constant mode C becomes dominated. Thus the key of obtaining the scale invariant spectrum 
during the radiation/matter domination is making the growing mode of $ spectrum have 
an opportunity to be inherited by the constant mode after the thermalisation. However 
unfortunately, it has been shown 3^ that in the simple scenario the growing mode of $ can 
hardly be matched to the constant model after the transition. This result is disappointing. 
However, it may be conceivable that our simplified operation in the calculation of primordial 
spectrum may have missed what. We here only phenomenally simulate the NEC violating 
process by using the scale field violating the NEC, and actually do not know the details 
of this NEC violating quantum fluctuation and subsequent thermalisation. The latter may 
significantly affect the final matching result. The example that in the gauge, in which 
instead of occurring a constant energy hypersurface the thermalisation may be everywhere 
simultaneous, the growing mode of $ spectrum can be inherit by the constant mode at late 
time has been pointed out j3, Q 3 . We will tentatively take this optimistic matching for 
the following calculations and leave behind some further comments in the final conclusion. 
Because in superhorizon scale the Bardeen potential is increased during the NEC violating 
fluctuation and constant after the end of fluctuation, it may be reasonable to take the value 
at the end of fluctuation to calculate the amplitude of perturbation, 

k 3/2 $ k - - Vy 2 /2- (18) 

where (JTTJ) has been used. 

The observations of CMB constrain the amplitude of perturbation is /c 3//2 Cfc ~ k 3 ^ 2 ^^ ~ 
10~ 5 , which require V e ~ 10~ 10 , which corresponds to h 2 e ~ 10~ 10 /|e|. We assume the 
thermalisation is almost instantaneous, thus h 2 e ~ p r , where p r is the energy density of 
radiation after the thermalisation. For example, for p r ~ m^ w ~ 10~ 60 , where m ew is the 
electroweak scale, which may be regarded as a loose lower limit of reheating temperature 
though the lowest limit might be in nucleosynthesis scale, we obtain |e| ~ 10 50 . The change of 
a is far smaller than that of h, thus from (j3J), we have J\f ~ In (h e /hi) ~ In (Vl^/jel^A 1 / 2 ) ~ 
127 — In |e|/2 ~ 69, which is enough to solve the horizon problem of our observable universe 



3 The recent numerical studies jsfii Ist| shown that the comoving curvature perturbation £ passes continuous 
through the transition and therefore the spectrum inherited by the Bardeen potential at late times is the 
one of C and not the spectrum of the growing mode of the Bardeen potential. However, the examples (or 
classes) discussed in these studies are still limited, thus whether the result is universal remains open. 
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42j . These results are also implicitly reflected in the depiction of causal structure of island 
universe of Fig. 2. . 

How do we distinguish the island universe from the inflation model by the observations? 
For e — > — oo, from (|13j) . we have v ~ 1/2, and thus the index of tensor perturbation is 2, 
which means that the island universe model will produce a more blue gravitational wave 
spectrum than the inflation model. This character is also actually a reflection of rapid change 
of background during the NEC violating fluctuation, compared with the nearly unchanged 
background during the inflation. This result generally leads to the intense suppression of 
gravitational wave amplitude on large scale. Thus it seems for the island universe model 
to be hardly possible to search for the imprint of gravitational wave in the CMB, while a 
stochastic gravitational wave will be consistent with the inflation model. 

In conclusion, it seems to be possible that we can live in "Islands", but slightly not 
optimistic, since the loophole how the perturbations propagate through the thermalisation 
surface remains, which may be fatal to the island universe model. However, it should be 
mentioned that the NEC violating fluctuation is actually a quantum phenomenon, and 
thus the innocence about the thermalisation and matching hypersurface may to some extent 
reflect the lack of our understanding on the quantum characters of gravity. We think that (at 
least) at present it will be premature to either confirm or deny the island universe model. 
The deep insights into the quantum gravity in the future may show us the final answer. 
However, whatever the outcome, it appears that the viability of the island universe model 
hangs in the balance and the different outcomes will either hold out or end its competing 
as an alternative scenario of the origin of observable universe. We will back to the relevant 
issues in the future. 
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